FOURTH-ORDER BESSEL-TYPE SPECIAL FUNCTIONS: 

A SURVEY 



W.N. EVERITT 

This paper is dedicated to the memory and achievements of 
George Neville Watson (1886 to 1965) 



Abstract. This survey paper reports on the properties of the fourth-order Bessel-type 
linear ordinary differential equation, on the generated self-adjoint differential operators in 
two associated Hilbert function spaces, and on the generalisation of the classical Hankcl 
integral transform. 

These results are based upon the properties of the classical Bessel and Laguerre second- 
order differential equations, and on the fourth-order Laguerre-type differential equation. 
From these differential equations and their solutions, limit processes yield the fourth-order 
Bessel-type functions and the associated differential equation. 
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1. Introduction 

This survey paper is based on joint work with the following named colleagues: 

Jyoti Das, University of Calcutta, India 
D.B. Hinton, University of Tennessee, USA 

H. Kalf, University of Munich, Germany 
L.L. Littlejohn, Utah State University, USA 
C. Markett, Technical University of Aachen, Germany 
M. Plum, University of Karlsruhe, Germany 
M. van Hoeij, Florida State University, USA 



2. History 

We see below that the structured definition of the general-even order Bessel-type special 
functions is dependent upon the Jacobi and Laguerre classical orthogonal polynomials, and 
the Jacobi-type and Laguerre-type orthogonal polynomials. 

These latter orthogonal polynomials were first defined by H.L. Krall in 1940, see [THl and 
|14j . and later studied in detail by A.M. Krall in 1981, see and by Koornwinder in 1984, 
see |[0|- In this respect see the two survey papers, [7j of 1990 and |Hj of 1999. 

The Bessel-type special functions of general even-order were introduced by Everitt and 
Markett in 1994, see 0. 

The properties of the fourth-order Bessel-type functions have been studied by the present 
author and the seven colleagues named in Section [T] above, in the papers j2], jlj and j^j. 



3. The fourth-order differential equation 

The fourth-order Bessel-type differential equation takes the form 

(3.1) (xy"(x))" - ((92T 1 + 8M- 1 x)y'{x)) / = Axy(x) for all x G (0, oo) 

where M G (0, oo) is a positive parameter and A G C, the complex field, is a spectral 
parameter. The differential equation (|3.1|) is derived in the paper jSJ Section 1, (1.10a)], by 
Everitt and Markett. 

This linear, ordinary differential equation on the interval (0, oo) c R, the real field, is 
written in Lagrange symmetric (formally self-adjoint) form, or equivalently Naimark form, 
see [HI Chapter V]. 

The structured Bessel-type functions of all even orders, and their associated linear dif- 
ferential equations, were introduced in the paper jSJ Section 1] through linear combinations 
of, and limit processes applied to, the Laguerre and Laguerre-type orthogonal polynomials, 
and to the classical Bessel functions. This process is best illustrated through the following 
diagram, see (HI Section 1, Page 328] (for the first two lines of this table see the earlier work 
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of Koornwinder and Markett JTj ) : 



(3.2) 



Jacobi polynomials 

k(a,(3)(l - x) a {l + xf 

i 

Laguerre polynomials 
k(a)x a exp(— x) 
I 

Bessel functions 

n(a)x 2a+1 



Jacobi-type polynomials 

k(a, f3)(l - x) a (l + x) p 
+MS{x + 1) + N5{x - 1) 
i 

Laguerre-type polynomials 
k(a)x a exp(— x) + NS(x) 
i 

Bessel-type functions 



K[a)x 



2a+l 



+ M8(x) 



The symbol entry (here k and k are positive numbers depending only on the parameters a 
and (3) under each special function indicates a non-negative (generalised) "weight", on the 
interval (—1, 1) or (0, oo), involved in: 

(a) the orthogonality property of the special functions 

(6) the weight coefficient in the associated differential equations. 

It is important to note in this diagram that: 

(z) a horizontal arrow — > indicates a definition process either by a linear combination 

of special functions of the same type but of different orders, or by a linear-differential 

combination of special functions of the same type and order (alternatively by an 

application of the Darboux transform, see |l(Jj ) 
(ii) a vertical arrow J, indicates a confluent limit process of one special function to give 

another special function 
(Hi) the use of the symbol MS(-) is a notational device to indicate that the monotonic 

function on the real line IR defining the weight has a jump at an end-point of the 

interval concerned, of magnitude M > 
(iv) the combination of any vertical arrow j with a horizontal arrow — > must give a 

consistent single entry. 

Information about the Jacobi-type and Laguerre-type orthogonal polynomials, and their 
associated differential equations, is given in the Everitt and Littlejohn survey paper [7]; see 
in particular the references in this paper to the introduction of the fourth-order Laguerre- 
type differential equation by H.L. and A.M. Krall, Koornwinder and by Littlejohn. The 
general Laguerre-type differential equation is introduced in the paper ^T] by Koekoek and 
Koekoek; the order of this linear differential equation is determined by 4 + 2a with a G No = 
{0,1,2,-..}. 

It is significant that the general order Bessel-type functions also satisfy a linear differential 
equation of order 4 + 2a (with a 6 No), being an inheritance from the order of the general 
Laguerre-type equation. 

The purpose of this survey paper is to discuss the properties of the Bessel-type linear 
differential equation in the special case when a = 0, as given in the bottom right-hand 
corner of the diagram; this is the fourth-order differential equation (J3.1|) and involves the 
weight coefficient n(0)x; its solutions should, in some sense, have orthogonality properties 
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with respect to the generalised weight function n(Q)x + M5(Q), where M > is the parameter 
appearing in the differential equation (|3.1J) : see [El Section 4]. 

Our knowledge of the special function solutions of the Bessel-type differential equation 
(J3.1j) is now more complete than at the time the paper [8] was written. However, the results 
in [HI Section 1, (1.8a)], with a = 0, show that the function defined by 

(3.3) J°x M {x) := [1 + M(A/2) 2 ] J {Xx) - 2M{X/2) 2 {Xx)- 1 J ± (Xx) for all x G (0, oo), 

is a solution of the differential equation (j3.1|) . for all A G C, and hence for all A G C, and all 
M > 0. Here: 

(i) the parameter M > 
(ii) the parameter A G C 

(Hi) the spectral parameter A and the parameter M, in the equation (J3.1j) . and the pa- 
rameters M and A, in the definition (J3.3j) . are connected by the relationship 

(3.4) A = A(A, M) = A 2 (A 2 + 8M" 1 ) for all A G C and all M > 

(iv) Jo and J\ are the classical Bessel functions (of the first kind), see [2H Chapter III]. 
Similar arguments to the methods given in [S] show that the function defined by 

(3.5) Y^ M (x) : = [1 + M(A/2) 2 ]F (Ax) - 2M(A/2) 2 (Ax)- 1 F 1 (Ax) for all x G (0, oo), 

is also a solution of the differential equation (J3.1|) . for all A G C, and hence for all A G C and 
all M > 0; here, again, Y and Y\ are classical Bessel functions (of the second kind), see [211 
Chapter III]. 

The earlier studies of the fourth-order differential equation (|3.1|) failed to find any ex- 
plicit form of two linearly independent solutions, additional to the solutions J^ M and Y®' M . 
However, results of van Hoeij, see [22] and [2B], using the computer algebra program Maple 
have yielded the required two additional solutions, here given the notation I®' M and K^ M , 
with explicit representation in terms of the classical modified Bessel functions Jo, Kq and 
h,Ki. These two additional solutions are defined as follows, where as far as possible we 
have followed the notation used for the solutions J®' M and Y®' M , 

(i) given A G C, with arg(A) G [0, 2ir), M G (0, oo) and using the principal value of yf, 
define 

(3.6) c = c(A, M) := V\ 2 + 8M- 1 and d = d(X, M) := 1 + M(A/2) 2 
(ii) define the solution /°' M , for all x G (0, oo), 

I° x ' M (x) := -dl (cx) + IcMx^I ^cx) 

(3.7) X := -[l + M(A/2) 2 ]Jo(xV A 2 + 8M - 1 ) 

+ y/M(2 + M(X/2) 2 )x- 1 h (xVX 2 + 8M- 1 ) 

(in) define the solution K®' M , for all x G (0, oo), 

I° x ' M (x) := dK (cx) + \cMx- x K x (cx) 

(3.8) ' := [1 + M(X/2) 2 }K (xV X 2 + 8M' 1 ) 

+ yjM(2 + M(X/2) 2 )x- 1 K 1 (xVA 2 + 8M" 1 ) 

Remark 3.1. We have 
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(1) The four linearly independent solutions J°' M , Y® , , i^' M provide a basis for all 
solutions of the original differential equation 1)3. subject to the (A, A) connection 
given in (J3.4j) . 

(2) These four solutions are real-valued on their domain (0, oo) for all A G R. 

(3) The domain (0, oo) of the solutions J°' and I x ' can be extended to the closed 
half-line [0, oo) with the properties 

J°' M (0) = J°' M (0) = 1 for all A G R and all M G (0, oo). 

The classical Bessel differential equation, with order a = 0, written in a form comparable 
to the fourth-order equation 1)3.10 . is best taken from the left-hand bottom corner of the 
diagram (J3.2|) ; from [HI Section 1, (1.2)] with a = we obtain 

(3.9) -(xy'(x))' = X 2 xy{x) for all Z G (0, oo); 

here A G C is the spectral parameter. It is to be observed that, formally, if the fourth-order 
Bessel-type equation (j3.1|) is multiplied by the parameter M > and then M tends to zero, 
we obtain essentially the classical Bessel equation of order zero (J3.9|) . on using the spectral 
relationship (J3.4)) between the parameters A and A. This Bessel differential equation ()3.9|) 
has solutions Ji(Ax) and Yi(Xx) for all x G (0, oo) and all A G C. 

For the need to apply the Frobenius series method of solution we also consider the differ- 
ential equation (J3.1)) on the complex plane C: 



(3.10) 



w^{z) + 2z- 1 w^\z) - (9z~ 2 + %M- r )w"(z) 
+(9z- 3 - %M- 1 z- l )w'{z) - Aw{z) = 



for all z G C. In this form the equation has a regular singularity at the origin 0, and an 
irregular singularity at the point at infinity oo of the complex plane C; all other points of 
the plane are regular or ordinary points for the differential equation. It should be noted that 
the classical Bessel differential equation (J3.9)) has the same classification when considered in 
the complex plane C. 

A calculation shows that the Frobenius indicial roots for the regular singularity of the 
differential equation ()3.10|) at the origin 0, are {4, 2, 0, —2}. The application of the Frobenius 
series method, using the computer programs [T| and Maple (see [22]), yield four linearly 
independent series solutions of ()3.1Uj) . each with infinite radius of convergence in the complex 
plane C. If these solutions are labelled to hold for the Bessel-type differential equation ()3.1|) 
then we have four solutions {y r {-, A, M) : r = 4, 2, 0, —2}, to accord with the indicial roots, 
to give the theorem, see [3 Section 3]: 

Theorem 3.1. For all A G C and all M > 0, the differential equation (J J. i|) has four linearly 
independent solutions {y r {-, A, M) : r = 4, 2, 0, —2}, defined on (0, oo) x C, with the following 
series properties as x — ► + , where the O -terms depend upon the complex spectral parameter 
A and the parameter M, 

{y±(x, A, M) = x 4 + lM~ 1 x 6 + 0(x 8 ) 
y 2 (x,A,M) = kx 2 + 0(x 4 \\n(x)\) 
y (x,A,M) = Z + 0(x 4 |ln(x)|) 
y- 2 {x,A,M) = rax' 2 + 0(\ln(x)\). 
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Here the fixed numbers k,l,m G M and are independent of the parameters A and M; these 
numbers are produced by the Frobenius computer program and have the explicit values: 

k = -(27720)" 1 , I = (174636000) -1 , m = -(9779616000)" 1 . 

4. Higher-order differential equations 

As mentioned in Section 01 above there exist Bessel-type linear differential equations of 
all even-orders 4 + 2a, where a G No is any non-negative integer. The definition and 
some properties of these differential equations, and the associated Bessel-type functions, are 
considered in detail in Sections 2 and 3]. 

Here we give the form of the sixth-order and eighth-order differential equations, as given in 
[HI Section 1, (1.10b) and (1.10c).]. (Note that there is a printing error in the display (1.10b); 
the numerical factor 255 is to be replaced by 225. Also printing errors in the display (1.10c) 
which are now to be corrected using the form of the differential equation ()4.2|) below.) 

(i) The sixth-order equation derived from the corrected differential expression for 
Section 1, (1.9) and (1.10b)] is 

, v -(x 3 y( 3 \x))W + (33xy"(z))" - ((2250T 1 - 96M~ 1 x 3 )y'(x))' 

[ ' = (A 6 + M~ l 2 A (3!) X 2 )x 3 y(x) for all x G (0, oo), 

where, as before, the parameters M G (0, oo) and A G C. 

When this equation is considered in the complex plane C the Frobenius indicial 
roots for the regular singularity at the origin are {6,4,2,0,-2,-4}, using the 
methods provided by 

(ii) The eighth-order equation derived from the corrected differential expression for |HJ 
Section 1, (1.9) and (1.10c)] is 

(x 5 y (4 \x))^ - (78x 3 y^(x))^ + (1809xy" (x))" 

(4.2) -((11025a;- 1 - 2 6 (4!)M-Vy(x))' 

= (A 8 + M" 1 2 6 (4!)A 2 )a; 5 ?/(x) for all x G (0, oo), 

where, as before, the parameters M G (0, oo) and A G C. 

When this equation is considered in the complex plane C the Frobenius indicial 
roots for the regular singularity at the origin are {8,6,4,2,0,-2,-4,-6}, using 
the methods provided by |2^] . 

We note that, formally, if the equations in (j4.1|) and (|4.2|) are multiplied by M > 0, and 
then letting M tend to zero we obtain, respectively, the two Sturm-Liouville differential 
equations, see [EI Section 1, (1.2)], 

(4.3) -(x 3 y'(x))' = \ 2 x 3 y(x) and - (x 5 y'(x))' = \ 2 x 5 y(x) for all x G (0, oo). 

For the solutions of these equations in classical Bessel functions see Section 1, (1.4)]. 

5. The fourth-order differential expression L m 
We define the differential expression Lm with domain D(Lm) as follows: 
(5.1) D(L M ) := {/ : (0, oo) - C : /W G AC loc (0, oo) for r = 0, 1, 2, 3}, 
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and for all / G D(Lm) 

(5.2) L M [f](x) := (xf'(x))" - {(Qx- 1 + 8M- 1 x)f(x))' (x G (O.oo)); 
it follows that 

(5.3) L M :D(L M )^Ll c (0,oo). 

The Green's formula for Lm on any compact interval [a,0\ C (0, +oo) is given by 

(5-4) I {g(x)L M [f](x)-f(x)l^g\(x)}dx = [f,g](x)f a , 



where the symplectic form [•, •](■) : D(Lm) x D(Lm) x (0, +oo) — > C is defined by 
\f,9](?) ■= g(x)(xf"(x))' - (xg"(x))'f(x) 
-x(g'(x)f"(x)-g"(x)f'(x)) 
(5.5) - (92T 1 + 8M _1 a;) (g(x)f'(x) - g'(x)f(x)) . 

The Dirichlet formula for L M on any compact interval [a, j3] C (0, +oo) is given by 

{xf"(x)g"(x) + (9X' 1 + 8M~ l x) f'(x)g'(x)}dx 



(5.6) 



I/ 3 



[f,g] D (x)f a + / L M [f](x)g(x)dx, 



where the Dirichlet form [-,-]d '■ D(Lm) X Dq(Lm) x (0, +oo) — > C is defined by, for / G 
D(L M ) and # G D (L M ), with 

(5.7) ^o(Im) := : (0, +oo) -> C : g {r) G Al7 1oc (0, +oo) for r = 0, 1} 
and 

(5.8) [/, g] D (x) := -^(x) + g'{x)xf"{x) + g(x) {9x^ + 8M~ l x) f'(x). 

6. HlLBERT FUNCTION SPACES 

The spectral properties of the fourth-order Bessel differential equation 

(6.1) {xy"{x))" - ((90T 1 + 8M~ 1 x)y'(x))' = Axy(x) for all x G (0, oo), 

with A G C as the spectral parameter, are considered in two Hilbert function spaces: 
(1) The Lebesgue weighted space 

(6.2) L 2 ((0,oo);x) := j/ : (0, +oo) -> C : jf x \f(x)\ 2 dx < +oo| 
with inner-product and norm defined by, for all /, g G £ 2 ((0, oo); x), 

/■oo 

(6.3) (/,<?):=/ dx and ||/|| := (/, J) 1 ^. 

Jo 

This space takes into account the weight function x on the right-hand side of (jfi.lj) . 
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(2) The Lebesgue-Stieltjes jump space L 2 ([0, oo); rrik), as suggested by the results in jHJ 
Section 4]. 

Let the monotonic non- decreasing function rhk : [0, oo) — > [0, oo) be defined by, 
where k > is a real parameter, 

"rhk{x) = —k for x = 

= x 2 /2 for all x G (0, +oo). 

Then generates a Baire measure on the a-algebra i3 of Borel sets on the 
interval [0, oo); in turn this measure generates a Lebesgue-Stieltjes integral for Borel 
measurable functions. 

The Hilbert function space L 2 ([0, oo); rrik) is defined on all functions with the prop- 
erties: 

(i) f : [0, oo) — > C and is Borel measurable on [0, oo) 

f°° 2 

(ii) / x\f(x)\ dx < +oo. 
Jo 

The norm and inner-product in L 2 ([0, oo); rrik) are defined by 

p roc 

(6.4) 11/11?:=/ \f(x)\ 2 dm k (x) = k\f(0)\ 2 + / x\f(x)\ 2 dx 

J[0,oa) JO 

and 

(6.5) (f,g)k-= f(x)g(x) dm k (x) = kf(0)g(0)+ xf(x)g(x) dx. 

J[0,cx>) JO 

Note that the first integrals in both these definitions are Lebesgue-Stieltjes integrals 
taken over the set [0, oo), whilst the second integrals can be taken as Lebesgue inte- 
grals. 



7. Differential operators generated by L 



M 



The Lagrange symmetric differential expression Lm generates self-adjoint operators in 
both the Hilbert function spaces £ 2 ((0, oo); x), and in L 2 ([0, oo); rrik) for all k G (0, oo). 

In the space L 2 ((0, oo); x) the expression Lm generates a continuum {T} of self-adjoint 
operators, including the significant Friedrichs operator F; these properties are developed and 
considered in Sections IH1 to [T^l below. 

For each k G (0, oo) the expression Lm generates a unique self-adjoint operator S k in the 
space L 2 ([0, oo);m fc ); the properties of this operator are considered in Sections IT41 and Hoi 

8. Differential operators in L 2 ((0,oo);x) 

The maximal and the minimal differential operators, denoted respectively T\ and T , as 
generated by the differential expression Lm in the Hilbert function space L 2 ((0, oo); x), are 
defined as follows, see [TH1 Chapter V, Section 17]: 

(i) T\ : D(Ti) C L 2 ((0, oo); x) - L 2 ((0, oo); x) by 
(8.1) D(Ti) := {f G D(L M ) : /, x^LmU) e ^((0, oo); x)} 
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and 

(8.2) TJ := x^LmU) for all / G D(T X ). 
From the Green's formula ()5.4|) if follows that the limits 

(8.3) [f,g](0 + ) := \ira[f,g](x) and [/, <?](oc) := lim [f,g](x) 

both exist and are finite in C for all f,gE D(T{). 
(ii) T : D{T ) C T 2 ((0, oo); x) -> L 2 ((0, oo); x) by 

D(T ) := {/ G D(Tx) : lim^of/, <7](z) = and 
K } lim^oof/^Kar) = for all /,<? G D(Tx)}, 

and 

(8.5) T / := rc- 1 L M (/) for all / G D(T ). 

From standard results we have the operator properties, see [THJ Chapter V], 

(8.6) To C Ti, T * = T x and T* = T , 

thereby noting that both T and T x are closed linear operators in T 2 ((0, oo); x). 

9. Self- adjoint operators in T 2 ((0, oo);x) 

In the weighted space T 2 ((0, oo); x) the Lagrange symmetric (formally self-adjoint) differ- 
ential expression has the following endpoint classifications at the singular endpoints and 
+oo (for additional details see Section 6]): 

(i) At + the singular endpoint is limit-3 in T 2 ((0, oo);x) 

(ii) At +oo the singular endpoint is Dirichlet and strong limit-2 in T 2 ((0, oo); x). 

Based on this information the self-adjoint extensions of the closed symmetric operator 
To are determined by the GKN-theorem on singular boundary conditions as given in ^HJ 
Chapter V] and [HJ. In particular, for the operators T and T 1; any self-adjoint operator T = 
T* generated by L M in T 2 ((0, oo);x) is a one-dimensional extension of T or, equivalently, a 
one-dimensonal restriction of T\. Let the domain D(T) as a restriction of the domain D(Ti) 
be determined by 

(9.1) D(T):={feD(T 1 ):[f, V }(0)=0} 1 

where the function ip G D(T\) is a non-null element of the quotient space D(Tx)/D(Tq) 
which satisfies the GKN symmetry condition 

(9.2) b,^](0) = 0. 
Then the differential operator T defined by 

(9.3) Tf := x- l L M [f] for all / G D{T) 

satisfies T* = T, and is self-adjoint in the Hilbert space T 2 ((0, oo); x). All such self-adjoint 
operators are determined in this way on making an appropriate choice of the boundary 
condition function ip. 
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10. Boundary properties at + 

The results of the following theorem are essential to obtaining the explicit forms of the 
boundary conditions at + to determine all self-adjoint extensions of To. 

Theorem 10.1. Let f E D(T\); then the values of /,/',/" can be defined at the point so 
that the following results hold: 

(i) f e AC[0, 1] 

(it) f E AC[0, 1] and f (0) = 

(iii) f" EAC loc (0,l] andf"eC[0,l] 

(iv) f^ e AC loc (0, 1] and \im x ^ 0+ (xf^(x)) = 0. 

For the proof of this theorem see [2J Section 8]. 

We consider the functions l,x, x 2 on the interval [0,1] but "patched", see the Naimark 
patching lemma ^Hl Chapter V, Section 17.3, Lemma 2], to zero on [2, oo) in such a manner 
that the patched functions belong to the domain D(Ti); we continue to use the symbols 
1, x, x 2 for the patched functions. 

A calculation shows that the results given in the next lemma are satisfied: 

Lemma 10.1. The patched functions l,x, x 2 have the following limit properties in respect 
of the symplectic form [•, •] and the maximal domain D(T{): 

(i) l,x 2 E D(T X ) butx£ D{T X ) 
(it) [1, 1](0+) = [x,x}(0+) = [x 2 ,x 2 ](0+) = 

(iii) [x,x 2 }(0 + ) = and [l,x 2 ](0 + ) = 16 

(iv) [l,x](0 + ) does not exist. 

The lemmas and corollaries now given below are taken from [21 Section 9], where proofs 
are given in detail. 

The results of Theorem 1 1 . 1 1 and Lemma HO.ll now provide a basis for the two-dimensional 
quotient space D(Ti)/D(T ); 

(10.1) D(Ti)/D(T ) = span{l, x 2 } = {a + bx 2 : a,b E C}. 

The linear independence of the functions {l,x 2 } within the the quotient space follows from 
the property [l,x 2 ](0 + ) = 16 ^ 0. 

A calculation now gives, recall Theorem ll0.ll 

Lemma 10.2. Let f £ D(Tx); then the following identities hold: 

W [/,l](0 + ) = -8f'(0) 

(ii) [f,x 2 }(0+) = 16/(0). 

Similarly we have 
Lemma 10.3. Let f,g e D(T\)\ then 

(i) lf,g](0 + ) = 8[f(0)g"(0)-f"(0)g(0)] 
(H) [f,g] D (0 + )=8f"(0)g(0)- 

We have the corollaries: 
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Corollary 10.1. The domain of the minimal operator Tq is determined explicitly by 

(10.2) D(T ) = {fe D{T X ) : /(0) = and /"(0) = 0}. 
Corollary 10.2. For all f G D(T X ) 

(10.3) J {x\f"(x)\ 2 + (9X- 1 + 8M- 1 x)\f'(x)\ 2 }dx <oo. 
Corollary 10.3. For all f,gE D(T X ) the Dirichlet formula takes the form 

POO 

(10.4) [T x f,g) = 8f (0)^(0) + / {xf"(x)g"(x) + (93T 1 + SM^x) f'{x)g\x)} dx. 

Jo 

11. Explicit boundary condition functions at + 

We can now determine all forms of the boundary condition function ip satisfying the 
symmetry condition (|9.2j) to determine the domain of all self-adjoint extensions T of the 
minimal operator T . 

Lemma 11.1. All self-adjoint extensions T of T generated by the differential expression 
Lm in L 2 ((0,oo;x) are determined by, using the patched functions l,x 2 , 

(11.1) D(T) := {/ G D(Ti) : [/, y?](0 + ) = where (i) (p(x) = a + (3x 2 

(ii) a, f3 eR and a 2 + (3 2 y^0}. 

and 

(11.2) (Tf) (x) := x' l L M (f)(x) for all x G (0, oo) and all f G D{T). 

There is an equivalent form of this last result, using the results of Lemma flO. 21 

Lemma 11.2. All self-adjoint extensions T of T generated by the differential expression 
Lm in L 2 ((0,oo;x) are determined by 

(11.3) D(T) := {/ G D(T{) : (z) -a/"(0) + 2/3/(0) = 

(ii) a,j9e R anrfa 2 + /3V 0}. 

and 

(11.4) (T/) (x) := x- l L M {f){x) for all x G (0, oo) and all f G D(T). 

Remark 11.1. We note the two special cases: 

(i) When a = the boundary condition is /(0) = 0; this boundary condition plays a 
special role, and gives an explicit form of the domain of the Friedrichs extension F 
of T ; see Section ITBI below. 

(ii) When (3 = the boundary condition is f"(0) = 0. 
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12. Spectral properties of the fourth-order Bessel-type operators 

Theorem 12.1. The minimal operator T , defined in and is bounded below in 

the space L 2 ((0, oo); x) by the null operator O, i.e. 

(12.1) (T f,f)>OforallfeD(T ). 

Proof. Since T is a restriction of the maximal operator T\ the result of Corollary 111). 31 can 
be applied to give, using also Corollary 110. 1\ 

POD 

(T /, /) = 8f (0)7(0) + / {x \f"(x)\ 2 + (9X- 1 + 8M-*x) \f'(x)\ 2 } dx 

Jo 

{x \f"(x)\ 2 + (90T 1 + SM^x) \f'(x)\ 2 } dx > 

for all f eD(T ). □ 

Theorem 12.2. (1) Let T be a self-adjoint extension ofT ; then: 
(i) The essential spectrum a ess (T) is given by 

(12.2) a ess {T) = cr amt (T) = [0,oo). 

(ii) There are no embedded eigenvalues of T in the essential spectrum. 
(Hi) T has at most one eigenvalue; if this eigenvalue is present then it is simple and 
lies in the interval (— oo,0). 
(2) Every point G (— oo, 0) is the eigenvalue of some unique self-adjoint extension T of 
T . 

Proof. The proof of this theorem is given in detail in |2J Section 13]. □ 

13. The Friedrichs extension F 

The closed symmetric operator T is bounded below in L 2 ((0, oo); x), see Theorem 112.11 
and the general theory of such operators implies the existence of a distinguished self-adjoint 
extension F, called the Friedrichs extension of Tq. 
This Friedrichs operator has the properties: 
(z) T C F = F* C Ti 
(ii) D(F) = {fe D(Ti) : f(0) = 0} 
(Hi) The essential spectrum cr ess (T) is given by 

(13.1) a css (T) =<T cont (F) = [0,oo) 

(iv) F has no eigenvalues. 

For a discussion of the definition and properties of this Friedrichs extension see [2J Section 
15]. 

14. Self- adjoint operator S k in L 2 ([0, oo); m k ) 

In this section, given any k G (0, oo), we define the operator S k generated by the differential 
expression L M in the Hilbert function space L 2 ([0, oo); m^), where this space is defined in 
Section El above. 
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Definition 14.1. Let k G (0, oo) be given; then the operator S k 

(14.1) S k : D(S k ) C L 2 ([0, oo); m k ) -> L 2 ([0, oo); m fc ) 

is defined by {see and and Theorem MU '. 1\ for the definition and properties of the 

domain D(Ti) C L 2 ((0, oo); x)) 

(i) L>(£ fc ) := -D(Tx) 
(ii) /or a// / G D^) 

/ 142 n f (Skf)(x) := -8fc- 1 f'(0)/ora; = 

1 ' ' \ := x _1 L M [/](x) /or all x G (0, oo). 

Theorem 14.1. For all k G (0, oo): 

(i) The linear manifold D(Sk) is dense in L 2 ([0, oo); m^). 

(ii) The operator Sk is hermitian in L 2 ([0, oo); mk). 

(in) The operator Sk is symmetric in L 2 ([0, oo); m k ). 

(iv) The operator Sk is bounded below in L 2 ([0, oo); m k ) 

(14.3) (S k fJ)k>OforallfeD(S k ). 

For the proof of this theorem see Theorem 5.2]. 

Theorem 14.2. Let k G (0, oo) be given; then the symmetric operator Sk on the domain 
D(Sk) is self-adjoint in the Hilbert function space L 2 ([0, oo); m^). 

For the proof of this theorem see jU Theorem 5.4]. 

Theorem 14.3. Let k G (0, oo) be given; then the operator S k on the domain D(Sk) is 
the unique self-adjoint operator generated by the differential expression Lm in the Hilbert 
function space L 2 ([0, oo); m k ). 

For the proof of this theorem see j3J Theorem 5.5]. 

15. Spectral properties of the self-adjoint operator Sk 
The spectral properties of the self-adjoint operator Sk in L 2 ([0, oo); m k ) are given by 

Theorem 15.1. For any k G (0, oo) let the self-adjoint operator S k in L 2 ([0, oo); m k ) be 
defined as in Definition s^. 1\ above; then the spectrum o~(S k ) of Sk has the following properties: 

(i) Sk has no eigenvalues 
(ii) the essential spectrum of Sk is given by 

(15.1) cr ess (Sk) = cr cont (S k ) = [0, oo). 

For the proof of this theorem see j3J Theorem 6.1]. 
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16. Distributional orthogonality relationships 

Recall that from the properties of the classical Bessel function J we have the result that 
Jo(-) ^ J 2 ((0, oo); x). However from Section 1, (1.7)] we have the following distributional 
(Schwartzian) orthogonal relationship for the classical Bessel function Jo, in the space V of 
distributions, 

POO 

(16.1) A / xJ (Xx)J (fix) dx = 5(X — fi) for all A, /i G (0, oo); 

Jo 

here 5 G T>' is the Dirac delta distribution. This is the generalised orthogonality property 
for the solutions Jo of the classical Bessel differential equation, of order 0, given by (|3.9j) : 
this result mirrors the spectral properties of this equation, when considered on the half-line 
(0, oo), in the space J 2 ((0, oo); x); in particular the result that every self-adjoint extension 
T of the corresponding minimal operator T has the property a ess (T) = [0, oo). 

The distributional proof of (|16.1|) is discussed in the forthcoming paper [5], where the 
result is also related to the properties of infinite integrals of Bessel functions as originated 
by Hankel, see [22 Chapter XIII]. 

As above for the Bessel function Jo we have, from the explicit representation ()3.3|) . the 
fourth-order Bessel-type function J°' M ^ J 2 ([0, oo); m^) for all k,M G (0, oo). 

To obtain a distributional orthogonality for J°' M , given any M > 0, it is necessary to 
choose a special value of the parameter k, i.e. k = M/2. Then it is shown in jHJ Section 4, 
Corollary 4.3] that we have the following distributional (Schwartzian) orthogonal relationship 
for the fourth-order Bessel-type function J?' , in the space T>' of distributions, 



(16.2) 



A [1 + M(A/2)T 2 { I q xJ° x ' M (x)J^(x) dx + fMJ A °< M (0) JM'(0) 

= <5(A — /i) for all A, /i G (0, oo). 



The distributional proof of (J16.2)) is discussed in the forthcoming paper [5]. 

As a formal representation it follows that ()16.2|) may be written as, using the inner-product 
for the space L 2 ([0, oo); m M / 2 ) , 

(16.3) A [1 + M(A/2) 2 ]" 2 ( J \ M (-)> J i' M (-)) M2 = <KA-//) for all A,/i G (0,oo). 

As another connection between the classical Bessel ()3.9|) and the fourth-order Bessel-type 
()3.1|) differential equations it is to be noted that, formally, the orthogonality result ()16.2|) 
tends to the orthogonality result (jl6-H) . as the parameter M tends to zero. 



17. The generalised Hankel transform 

From the general theory of symmetric integrable-square transforms given in [201 Chapter 
VIII] one form of the classical Hankel transform, for the Bessel function J and working in 
the Hilbert function space L 2 ((0, oo; x), is: 
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(i) Let / G £ 2 ((0, oo); x) then the Hankel transform g G L 2 ((0, oo); s) is given by, for 
s G (0, oo), 

pOO 

;i7.1) <?(*) = / £J (*0/(0 ^ 

Jo 

with convergence of the integral in I/ 2 ((0, oo); s) 
(ii) With (7 G £ 2 ((0, 00); s) the inverse transform, to recover /, is given by, for x G (0, 00), 

17.2) f{x) = / sJ (a:s)</(s) ds 



with convergence of the integral in L 2 ((0, oo);x) 
(in) The Parseval relation holds between g and / 

pOO /*O0 

17.3) / x\f(x)\ 2 dx= / s |g(s)| 2 ds. 



There is also a direct convergence form of the Hankel transform which is best written as, 
starting with / G L 1 ((0, 00); x), 

poo poo 

(17.4) f(x) = / sJ (xs) ds / C J (sO/(O d£ 

Jo Jo 

with x G (0, 00). Here the integrals are Lebesgue or limits of Lebesgue integrals as discussed 
in [201 Chapter VIII]. 

There is an equivalent generalised Hankel transform involving the fourth-order Bessel-type 
function J°' M (-) and working now in the Hilbert function space L 2 ([0, 00); m M / 2 ) note again 
these results require the unique choice of k = M/2. 

The complete discussion of the following results for the generalised Hankel transform are 
to be found in the forthcoming paper jj]. 

To state these results the Lebesgue-Stieltjes Hilbert function space L 2 ((0, oo);n) is re- 
quired. Let the function h : [0, 00) — > [0, 00) be defined by 

(17.5) h(X) := \\ 2 [I + M(A/2) 2 ] _1 for all A G [0, 00); 
then 

n'(A) = A [1 + M(A/2) 2 ] " 2 > for all A G [0, oo) 

so that n is monotonic increasing on [0, oo) and generates a Baire measure on the cx-algebra 
E of Borel sets on the interval [0, oo). The Hilbert space L 2 ((0, oo); n) is then defined as the 
set of all Borel measurable complex-valued functions / on [0, oo) such that 

|/(A)| 2 dn(A)<+oo, 

[0,oo) 

with norm and inner-product defined by 

poo poo 

(17.6) \\ff n := \f(X)\ 2 dn(X)= |/(A)| 2 A [l + M(A/2) 2 ] " 2 dX 



(f,9)n= / f(X)g(X)dn(X) = / f(X)g(X)X [l + M(A/2) 2 ] dX. 

J [0,oo ) 
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Remark 17.1. This norm ||-|| and inner-product (-, •)„ for the space L 2 ((0, oo); n) are not 
to be confused with the norm ||-|| fc and inner-product (•,•)&> introduced in Sectional for the 
space L 2 ([0, oo); m^). 

We note that the weight function A i — > A [1 + M(A/2) 2 ] 2 in the integral in 1)17.6)1 is the 
factor in the distributional orthogonal relationships (J16.2)) and 1)16.3)1 . 

(1) The L 2 -theory of the generalised Hankel transform is given by the following results: 

Theorem 17.1. Let f E L 2 ([0, 00)5771^/2) • Then there exists exactly one function 
g E L 2 ((0, 00); n) with the property that 



17.7) 



'17.81 



\g{\)\ 2 aln{\) 



[0,00) 



\f{x)\ 2 dm M/2 (x) 



here g is defined by, for almost all A G (0, 00) 



(Fuf) (A) := g(X) 



Jx M (x)f(x)dm M/2 (x), 



[0,00) 



thereby defining also the generalised Hankel operator 
17.9) T M ■ L 2 ([0,oo);m M/2 ) -> L 2 ((0, 00); n). 

In addition g satisfies 



17.10) 



g{\)dn{\) = /(0). 



Remark 17.2. Note that the result ()17.8j) has to be interpreted as follows, in (i) 
and (ii): 



(ii) 



/ J\' M \ x )f( x )dmM/2(x) G L 2 ((0, oo);n) for all X G [0, 00) 

J[0,X] 



lim 



Jx M (x)f(x)dm M/2 (x) 



[0,X] 



dn(X) = 0. 



(Hi) Note that the Cauchy-Schwarz inequality shows that if g G L 2 ((0, oo);n) then 
g G ^((0,00);^. 

Theorem 17.2. Let g G L 2 ((0, 00); n). Then there is exactly one function f G 
L 2 ([0, 00); rriM/2) with the property that \17.1 ) is satisfied; here f is defined by 



17.111 



(Qug) (x) := f(x) := ^ 



g(\)dn(X) 



if x = 



J°' (x)g(X)dn(X) for x G (0, 00) 



^ JO 



thereby defining also the inverse generalised Hankel operator 
17.12) Q M : L 2 ((0, 00); n) -> L 2 ([0, 00); m M/2 ) . 

Remark 17.3. Note that the result fjl7.ll)) has to be interpreted as follows: 
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>") / J x M {x)g{\)dn(\) G L 2 ([0, oo); m M/2 ) for all A G (0, oo) 



I* 

(if) 

/(*) - / J A °' M (a;)^(A)^(A) 



lim 

A— >oo 











xdx = 0. 



(2) The direct convergence of the generalised Hankel transform is given by the following 
results: 

Theorem 17.3. Let 7 G (0, 00). If f : (0, 00) — > K. /ios £/ie property 
17.13) x 1 — ► v^/O) e ^(0, 00) 

and 25 0/ bounded variation in a neighbourhood 0/7, i/ien 



17.14) 



/OO / POO \ 

Ja ° ,M(7) U ) 



Let /x G (0, 00). 7/ p : (0, 00) — > K /ias t/ie property 

:i7.15) A __^_ eil(0 , M) 



(17.16) 



and 0/ bounded variation in a neighbourhood of //, i/iera 
±[<?(// + 0) + <7(/i-0)] 



f 00 

^°' M (^) ( / J A °' M (x)r/(A)^(A) ) dm M/2 (x). 



Remark 17.4. The integrals in Theorem 117.31 are either Lebesgue integrals or limits 
of such integrals over compact intervals of (0, 00). 

Corollary 17.1. (i) If 7 G (0, 00) is a point of continuity of the function f then 
(TTj , JTOl and ( [777! zmp/y 

/■oo 

{Qm^mMi) = / 7 A °' M ( 7 )(^M/)(A)dn(A) 
Jo 

/■oo 

= (M/2)/(0) / J A °' M (7)^(A)+ 



OO / /"OO 

0,M/ \ / / T 0,A/, 



^a' (7) ( / J\ (x)f(x)xdx J dn(A) 
(17-17) = /( 7 ), 



since 



poo 

;i7.18) / J° X ' M \v)dn{\) = for all r, G (0, 00). 

Jo 

(ii) 7/7 = then 

/"OO 

(M^m/))(0) = (M/2)/(0) / dn(A) = /(0) 
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since 

POO /POO 

J°' M (0) I J J° x ' M (x)f(x)xdx ) (///(A) 0. 



from Ity , the use of the Fubini integral theorem and noting that J°' M (0) = 1. 
(Hi) If fi G (0, oo) is a point of continuity of g then l \17.8\j , l \17.1 i| ) and ( i 7. i 6j ) imply 



(Fm(Gm9))(ij>) = g{n)- 

18. The Plum partial differential equation 

The Plum equation is a fourth-order linear partial differential in the Euclidean space R 2 
of two dimensions, derived from a linear partial differential expression which is connected 
with the fourth-order Bessel-type ordinary differential equation. 

If the Laplacian V 2 partial differential expression is written in polar co-ordinates 

(9 2 1 (9 1 d 2 

(18-1) V^ + ijUIiL 

then the Plum equation has the form, with u = u(r, 6), 

(18.2) V 4 m - 7V 2 m - %u = Au. 



r 



Here 7 > is determined by 7 = 8M _1 where M > is the parameter in the fourth-order 
Bessel equation (|3.1|) . and A G C is a spectral parameter. 

Written out the equation (|18.2|) becomes, see (THl Section 1, (2)], 

d 4 u 2d 3 u ( 1 \ d 2 u ( 1 7\ du 1 d 4 u 



+ --5-7 + — rTkr+h-:k + 



2 d 4 u 2 d 3 u ( 4 7 \ <9 2 m 4 7 ^ 

(18.3) = Am. 

From the results given in JT!5] and jinj assume that a solution for ()18.2|) is of the separated 
form 

(18.4) u(r,0) = v(r)w(8) 

where u> is required to be a solution of the second-order Sturm-Liouville differential equation 

(18.5) = 4w(0). 
Note that w is then of the general form 

(18.6) w(6) = A cos(2#) + B sin(2#) 

for, say, 9 G [0, tt] and scalars A, B. 

Also note that the factor 4 in the equation (|18.5|) is critical, and has to be fixed, for the 
separation method to be effective. 
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Substitution of (jl8.4j) into (j!8.3|) yields, see [TBJ Section 1] and [15) . 

v^\r) + V(r) + (-1 - 7 ) v"(r) + (i - v'(r) + ^v(rf) w{9) 

+ (~- 2 v"(r) + I„'(r) + f-H + ^) «(r) - %(r)) w{9) 

(18.7) = At>(r)w(#) for all r G (0, oo) and all 9 G [0, tt]. 

For (jl8.7|) to hold requires that the function v (■), on gathering up terms, has to satisfy 
the ordinary differential equation, see [TBI Section 1, (4)], for all r G (0, oo), 

(18.8) v^(r) + -v^{r) + ( — % - v"(r) + ( \ - -) v'(r) - Av(r) = 0. 



This last equation may be written in the Lagrange symmetric form 

(18.9) (rv"(r))" - ((9r _1 + ~/r)v'(r))' = Arv(r) for all r G (0, oo), 

which is the Bessel fourth-order differential equation (|3.1)) when 7 = 8M _1 . 

Thus separated solutions of the partial differential equation (|18.2|) can be written in the 
form 

(18.10) u(r, 9) = v(r)w(9) for all r G (0, 00) and 9 G [0, 2tt], 

where w(-) is any trigonometrical solution ()18.6|) of ()18.5|) . and v (■) is any solution of the 
fourth-order Bessel equation ()18.9|) for any choice of the spectral parameter A. 
Defining the partial differential expression P 7 [-], for 7 G (0, 00), by 

(18.11) P 7 [w] := ^V 4 m - 7V 2 m - 

it is shown in ^Hl that P 7 is a formally symmetric linear partial differential expression in 
L 2 (E 2 ), using polar co-ordinates (r,9). 

Some early studies indicate that there may be problems in applied mathematics, for which 
the partial differential equation Pj[u] = Au is involved in one or more of the associated 
mathematical models. 
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